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SUMMARY 

Mixing of coaxial streams is analyzed by application of integral techniques. An 
integrated turbulent kinetic energy (TKE) equation is solved simultaneously with the 
integral equations for the mean flow. Normalized TKE profile shapes are obtained from 
incompressible jet and shear layer experiments and are assumed to be applicable to all 
free turbulent flows. The shear stress at the midpoint of the mixing zone is assumed to 
be directly proportional to the local TKE, and dissipation is treated with a generalization 
of the model developed for isotropic turbulence. Although the analysis was developed 
for ducted flows, constant -pressure flows were approximated with the duct much larger 
than the jet. The axisymmetric flows under consideration have been predicted with 
reasonable accuracy. Fairly good results have also been obtained for the fully devel- 
oped two-dimensional shear layers, which were computed as thin layers at the boundary 
of a large circular jet. 


INTRODUCTION 

An extensive integral analysis of ducted turbulent mixing processes (fig. 1) has 
been developed at the Arnold Engineering Development Center (AEDC) (refs. 1 and 2). 

As usual with such analyses, the shape of the velocity profile is assumed, and the inte- 
gral form of the mean flow governing equations is used to compute the shear -layer 
growth rate and other dependent variables of the problem. In the integral analysis of 
references 1 and 2, the turbulent shear stress at the midpoint of the shear layer is 
computed by use of a model for the turbulent eddy viscosity. The integral analysis has 
been extended to flows with equilibrium chemical reactions, to flows which extend across 
the entire mixing duct, and to flows in which the inviscid portion of the jet flow must be 
treated with the method of characteristics. The method has also been applied to flows 
with embedded recirculation zones. The analytical framework has been developed for 
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quite complex flow situations, but a serious deficiency has always been the model used 
for the turbulent eddy viscosity, the Prandtl incompressible model with a compressibility 
correction similar to that proposed by Donaldson and Gray (ref. 3). This eddy -viscosity 
model is not adequate when the secondary stream velocity exceeds about 0.2 times the 
primary stream velocity, or when large pressure gradients exist in the flow field. Other 
eddy -viscosity models have been proposed which would perhaps yield better results for 
certain flows, but based on the extensive evaluation by Harsha (ref. 4), it is not likely that 
any eddy -viscosity model will be applicable to the entire range of free turbulent flows of 
interest. A fundamental problem with eddy -viscosity models is that the turbulent shear 
stress is related only to the local mean flow properties (this local dependence is true only 
in simple limiting flow situations); most free turbulent flows are characterized by signifi- 
cant history effects on the turbulent transport. 

Starting with the work of Bradshaw and associates (ref. 5) a few years ago, con- 
siderable effort has been devoted to development of the turbulent kinetic energy (TKE) 
methods for turbulent shear flows. In these methods, the turbulent shear stress is related 
to the kinetic energy of the turbulent motion, and the TKE governing equation is solved 
simultaneously with the mean flow governing equations; the TKE equation is, in effect, a 
governing equation for the turbulent shear stress. Two different methods have been used 
to relate the shear stress to the TKE in free shear flows. In the work at AEDC (refs. 4, 

6, and 7), Bradshaw's direct relationship between shear stress and TKE has been used. 
Other investigators (refs. 8 to 11) have related shear stress to TKE by defining an eddy 
viscosity which is the product of a length scale and the square root of the TKE. Although 
there appears to be experimental evidence (ref. 12) for the Bradshaw approach, both 
approaches are great improvements over earlier eddy -viscosity models in that turbulent 
shear flows are recognized to be indeed turbulent, that is, to consist of both mean and 
fluctuating components. The current TKE methods are undoubtedly oversimplified for 
the whole spectrum of turbulent flows in nature; however, Harsha’s work (ref. 4) has 
shown that the TKE approach is useful for a large class of shear flows which is commonly 
encountered in engineering applications. 

In the present study, the earlier integral approach for ducted flows (refs. 1 and 2) 
has been extended to include an integrated TKE equation. Because the TKE equation has 
been integrated across the entire shear layer, no model for the lateral diffusion of TKE 
needs to be specified. In addition, the relation between the TKE and the shear stress is 
specified only at the midpoint of the shear layer. These simplifications are achieved 
with a penalty — the shape of the lateral TKE profiles in the shear layer must be speci- 
fied. The TKE profile shapes have been obtained from incompressible experiments but 
have been used with reasonable success for flows with large density gradients. 
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The present integral method is limited to those flows in which the velocity profile 
shape is essentially fully developed (shape similar) throughout the flow field; that is, the 
initial boundary layers must be relatively thin. Deference must be made to the more 
powerful finite -difference TKE methods for those flows which have developing velocity 
profiles over a significant axial distance. In addition, the present integral method has 
been formulated for only axisymmetric flow, and the two-dimensional jet and wake flows 
have not been computed. However, the fully developed shear layers (test cases 1 to 3) 
have been computed as thin layers at the boundary of a very large axisymmetric jet. 

SYMBOLS 

aj constant in TKE shear -stress relation 

a 2 coefficient in dissipation term 

b mixing zone width 

c correction factor for a 2 

C mass fraction of elements from primary stream 

D diameter of primary stream nozzle 

H stagnation enthalpy 

k turbulent kinetic energy 

k normalized turbulent kinetic energy 

K exponent in velocity-concentration relation 

M Mach number 

M e external stream Mach number 

M 0 central stream Mach number 

n boundary-layer profile exponent 
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p 


static pressure 


r radial coordinate 

ri radius of inner mixing zone boundary 

r Q radius of primary stream nozzle 

r w duct wall radius 

R gas constant 

R-j turbulent Reynolds number 

T static temperature 

Tt stagnation temperature 

u axial velocity component 

u e outer stream velocity 

u Q initial primary stream velocity 

u'^ square of turbulent velocity component in x-direction 

v radial or transverse velocity 

v' 2 square of turbulent velocity component in r-direction 

w TKE profile parameter (eq. (11)) 

square of turbulent velocity component in circumferential direction 

W wake center-line velocity defect , -■ ~ Uc 

u e 

x axial coordinate 

x CO re length of first regime 
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y 


transverse coordinate 


r - r, 

mixing zone coordinate, - 

b 

a mass fraction of primary stream species on center line 

6 boundary -layer thickness 

e turbulent eddy viscosity 

p density 

r turbulent shear stress 

a spreading parameter for two-dimensional shear layer 

a 0 spreading parameter at a reference condition 

Subscripts: 

1 initial mixing station; high speed stream for two-dimensional shear layer 

2 low -speed stream for two-dimensional shear layer 

a inviscid secondary flow 

W boundary layer 

c center line 

f far field 

j inviscid primary flow 

m half -velocity control surface in mixing zone 

n near field 
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si 

shear layer 

w 

duct wall 


DEVELOPMENT OF ANALYSIS 

The first and second regimes of figure 1 will be considered in this paper. A region 
of inviscid secondary flow exists throughout the duct, and the turbulent mixing zone is 
free turbulent in nature; that is, the turbulent flow is not adjacent to the wall. The duct 
wall interacts with the turbulent shear layer only through the axial -pressure gradients 
which are imposed by its presence. 

Fundamental Assumptions 

The following principal assumptions have been used in developing the analysis: 

(1) The flow is axisymmetric. 

(2) All gases obey the perfect gas law. 

(3) The usual boundary -layer assumptions are used; that is, negligible radial 
pressure gradients, and so forth. 

(4) The inviscid portions of the primary and secondary flows are one -dimensional 
and isentropic. 

(5) The mixing layer is completely turbulent, and the initial boundary -layer thick- 
nesses at the initiation of mixing are very small compared with the length of the first 
regime (fig. 1). 

(6) The thickness of the nozzle lip separating the primary and secondary flows 
is negligible. 

(7) The viscous effects at the duct wall are negligible. 

(8) The normalized velocity profiles in the mixing layer are similar in shape at 
all axial stations and are represented by a cosine function. 

(9) The turbulent Prandtl and Lewis numbers are unity. 

(10) The TKE profile shapes in the shear layer which have been obtained from 
constant -density experiments are unaffected by density gradients in the shear layer. 

(11) The turbulent kinetic energy outside the shear layer boundaries is negligible. 
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Basic Integral Equations 


Nomenclature for the analysis is illustrated in figure 1. In this section, the nomen- 
clature is generally consistent with that of references 1 and 2; the results are described 
in terms of the recommended nomenclature for this conference. 


By integrating the boundary -layer differential equations, five basic integral equa- 
tions are obtained: (1) a continuity equation for the entire flow, (2) a momentum equation 

for the entire flow, (3) a momentum equation for the flow between the duct center line and 
a control surface arbitrarily located at the midpoint of the shear layer, (4) a jet species 
conservation equation for the entire flow, and (5) a turbulent kinetic energy equation for 
the shear layer. 

Continuity equation: 


j* W dr = "Pw v w r w = -Pa u a r w 

Overall momentum equation: 

r r w 9 / j r w^ d ^w 2 

i ^-( P u 2 )rdr=-_ -j--P a vV 


dr w 

dx 


dr 


)q fe 2 dx 

Half -radius momentum equation: 

• r™ „ . rr 

J 0 OA ' — •'O 

Jet species conservation equation: 

>r. 

J 0 


w 
dx 


r 4( pu2)r ar - “■» r 4 (pu>r dr = 


i w 4 <puc)r dr = ° 


T m r m 


r 2 dp 

x m F w 
2 dx 


( 1 ) 


( 2 ) 


(3) 


(4) 


The differential form of the TKE equation is 


pur — + pvr — = rr — + Diffusion - Dissipation 
9x 9r 3r 

where 

k = |(u’ 2 +v’ 2 + w’ 2 ) 


By integrating the TKE equation across the entire shear layer, the lateral diffusion term 
disappears. The dissipation is given by the usual relation developed for isotropic turbu- 
lence (ref. 4): 


Dissipation = 



b 
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Thus, the basic integral TKE equation is 


r b ar - r 


r^+b 


9u . 
r — r dr 
9r 


ag r*r i +b 

b J r . 


pk^^r dr 


(5) 


Velocity Profile 
The velocity profile is given by 

(6) 

I* “I* * 

where y and u c = u- in the first regime. The cosine profile is shown in fig- 

b J 

ure 2, along with the experimental results of Spencer (ref. 13) for a fully developed two- 
dimensional shear layer. At the control surface, r m (= r^ + the velocity is 

u m = |( u c + u a) 


Concentration and Enthalpy Profiles 


For unity Lewis number, the normalized total enthalpy and concentration profiles 
are essentially identical and are related to the velocity profile by 


C = 


H " H a 
h“Th^ 



(7) 


The parameter K has been introduced so that jet species can be conserved in 
variable pressure flows. For unity turbulent Prandtl number, equation (7) with K = 1 
is identically true for constant -pressure flow; the pressure gradients in all the flows con- 
sidered in this communication are negligible, and K remained unity for all computations. 
Therefore, equation (4) could have been deleted for these flows, since the assumption of 
K = 1 would have identically satisfied that equation. 


Turbulent Kinetic Energy Profiles 

As with the velocity profile, the shape of the TKE profiles must be specified. The 
near -field (first -regime) TKE profiles have been obtained from the experimental results 
of Spencer (ref. 13) and Liepmann and Laufer (ref. 14) for constant -density two- 
dimensional shear layers. The experimental results are shown in figure 3, along with 
the following analytical function which has been fitted to the data: 

k n =JS_= 0.51 - 0.51 cos JOL (0£y£0.45) (8a) 

11 k m 0.45 


584 


(0.45 = y = 1.0) (8b) 


kn = 0.51 - 0.51 cos[2ir(0.0909 + 0.9091y>) 


In fitting the data, the recent experiments of Spencer were given more weight than the 
Liepmann and Laufer experiments because of the improvement of hot-wire techniques in 
the past 20 years. 


In the far field, well downstream of the end of the potential core, the experiments 
of Pindell (ref. 15) and Wygnanski and Fiedler (ref. 16) on the constant -density axisym- 
metric jet into still air were used. The experimental TKE profiles are shown in fig- 
ure 4, along with the following analytical function which has been fitted to the data: 

kf = 0.5 + sin 1 + 1 y) (9) 

Obviously, the shape of the TKE profile cannot change dis continuously from k n to 
kf at the end of the first regime. It has been hypothesized that a transition region exists 
in which the TKE profile evolves from k n to kf. In this transition region, the TKE 
profile is assumed to be a linear combination of k n and kf 

k = wkf + (1 - w)k n (10) 


where w is an empirical function of x. The resulting family of TKE profiles is shown 
in figure 5. Because many features of jets scale with x/x core , it has been hypothesized 
that w can also be related to x/x core . The experiments of Pindell (ref. 15), Sami 
(refs. 17 and 18), and Bradshaw et al. (ref. 19) on the constant -density jet were used to 
develop the following empirical function for w: 


w = 1 - exp 


-1.09(-£ l\ 

\ core J 


( 11 ) 


The center -line TKE, normalized with the value at r m , is shown in figure 6 for the 
constant -density jet into still air. A value of x core = 4.66D was used in computing the 
curve shown in figure 6; this core length is predicted by the present theory for a jet with 
negligible initial boundary layer. The TKE profiles are seen to be essentially fully 


developed (w * 1) for x greater than 5x 


core* 


It should be noted that the value of x 


core 


is not prescribed for a particular flow but is a result of the first -regime solution. 


Even though the TKE profile shapes were obtained from relatively simple constant- 
density flows, the two-dimensional shear layer and the axisymmetric jet, it is assumed 
that these profile relationships apply to all flows. It should be noted that nothing has been 
stated about the level of TKE, but only that the TKE profiles, normalized with k m , are 
given by equations (8) to (11). 


Relation Between TKE and Shear Stress 

As the present integral theory has been formulated, the relationship between r 
and k must be specified only at the midpoint of the shear layer r m . This method is 
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distinctly different from finite -difference TKE methods, in which it is necessary to spec- 
ify the variation of t with k across the entire mixing zone. The linear relationship 

| T m | = a lP m k m (12) 


is used, with aj = 0.3, the value found by Harsha and Lee (ref. 11) in the high shear region 
of a variety of constant -density free turbulent flows. This value of a^ is also the same 
as that used by Bradshaw et al. (ref. 5) for turbulent boundary layers. In order to have 
the proper sign on r m , equation (12) is written as 

a l^m k m( u a “ u c) 


The Dissipation Parameter a 2 

In the extensive evaluation of his differential TKE method (ref. 4), Harsha used 
a 2 = 1.5. During the development of the present integral method, a 2 = 1.69 was found 
to give good results for constant -density two-dimensional shear layers and for the axisym- 
metric jet into still air. Other flows, however, were found to require significantly differ- 
ent values for a 2 if reasonably good mean flow predictions were to be achieved. Speci- 
fically, supersonic flows require a 2 to be larger than 1.69, and some flows with very 
high shear stress levels require a 2 to be less than 1.69. Finally, it was found that a 2 
could be correlated with the turbulent Reynolds number Rrp which is defined as 


Rrp = 


u c - u a b 

e m 


where e m is the local eddy viscosity at 


fied but is computed from 


€ m - 


'm 


n 8U 


Pm -5? 

m 


r m . Of course, the eddy viscosity is not speci- 


For the cosine profile, 


9u 

9r 


m 


77 ( u a ~ u c) 
2b 


and R-p may be written as 

„ *0 - u a f 

T ~ 2a l k m 


5.236 


( u c - u a)‘ 




Thus, it has been found that the dissipation coefficient a 2 can be related to the ratio of 
the mean flow velocity difference across the shear layer to the turbulent velocity fluctua- 
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tion level (characterized by k m ) in the layer. The a 2 function which has been devel- 
oped is shown in figure 7, along with the equations which describe the function. 

The approach used in the development of the a 2 - Rp function will be briefly 
described. For Bp > 145, the a 2 function was developed by computational experi- 
ments on the near field of unheated compressible air jets exhausting into still air. 
Warren's potential core length data (ref. 20) for experimental flows with thin initial 
boundary layers and 0.69 < Mj < 2.6 were used to establish approximately a few desired 
a 2 values. A direct correlation of a 2 with a characteristic Mach number for the shear 
layer was abandoned because such a correlation fails for two -stream supersonic flows 
such as the combustion flows reported in reference 1. The turbulent Reynolds number 
was finally found to correlate consistently the a 2 values for the preceding flows. 
Because Rp is related to the local turbulence characteristics in the shear layer, it is 
more appropriate as a dissipation parameter than some other parameter which is related 
only to the mean flow in the layer. When it is considered that only the near -field results 
for a few experimental flows were used in developing the a 2 function for Rp > 145, the 
overall performance of this part of the a 2 - R-p function (fig. 7) has been reasonably 
satisfactory for a variety of flows. The high Rp portion of the a 2 function is subject 
to further refinement, however, particularly for Rp > 300, that is, for fully developed 
single stream flows with Mj > 2.7. 

Experience with far -field predictions of jets in moving streams indicated that a 2 
should be somewhat less than 1.69 for Rp < 70. The function used for R T < 70 (fig. 7) 
is the first one tried, and no attempt has been made to improve it. 


The effects of density ratio caused by jet Mach number are adequately predicted by 
the a 2 - Rp function; however, prediction of the entire range of flows of interest is 
improved if small additional corrections to a 2 are made as a function of density ratio. 
Tentatively, the following corrections have been developed and used: 


a 2 = 


s( r t) 


where ag(Rp) is as shown in figure 7 and c is given by 

c = 0.984 + 0.016 2i 
p jl 


or 


c = 0.95 + 0.05 


R a T ta 

R j T .j 


( p al =• »jl) 


(Pal < Pjl) 


These density corrections to the basic ^ - Rp function are perhaps required because 
the empirical TKE profiles are inadequate for flows with large density gradients; this 
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point will not be resolved until detailed turbulence structure data are available for flows 
with large density gradients. 

The a 2 function as described has been used for all of the shear layer and jet 
computations presented in this paper. Experience has shown that wakes require some- 
what different dissipation than jet flows, and the axisymmetric wake computations (test 
cases 15 and 17) were made with a 2 = lo40. 


Turbulence Production 


The first term on the right-hand side of equation (5) represents the production of 
turbulence by the shear stress. For the boundary conditions of the shear layer, the pro- 
duction is equal but opposite in sign to the dissipation of mean flow mechanical energy: 



ir r i + 
2 J 0 


— (pu 3 )r dr 
9x N 


1 2 f ri+b a, x ^ dPw r r i +D , 

+ - u z \ — (pu)r dr \ ur dr 

2 a Jq Hy Jo 


dp 

dx J o 


r,+b 


(14) 


By substituting equation (14) into equation (5), the following form of the integral TKE 
equation is obtained: 

I,. 1 ^ puk)r dr * - \ j„ 1 dr + i u a 2 ] o r ‘ +b dr 



ur dr - 


*2 r r i +b 
b J ri 


pk'^r dr 


By using the TKE equation in this form, the shape of the shear -stress profile need not be 
specified. The turbulence production is related only to the dissipation of mean flow 
mechanical energy, which, in turn, is related to the mean flow profiles and the rate of 
growth of the shear layer. 


Solution Technique 

Sufficient information is available to transform equations (1) to (5) into a system of 
ordinary differential equations which is linear in the derivatives of the dependent vari- 
ables ^dp W/ /dx, etc.). This transformation procedure is described in detail in reference 2. 
After the system of linear equations is solved for the derivatives, the resulting five dif- 
ferential equations are then numerically integrated with a modified Euler technique (vari- 
able step size). An IBM 370/155 digital computer was used to obtain the numerical solu- 
tions; a typical flow -field solution required a computation time of approximately 2 minutes. 
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In the first regime, the dependent variables are p w , r^, b, K, and km. In the 
second regime, the dependent variables are p w , u c , b, K, and k m . 

Initial Conditions 

In order to integrate the system of differential equations, initial values must be 
specified for each of the dependent variables. The most critical of these initial values 
is that for k m . As usual with TKE methods, the convective terms cause the initial 
condition for k m (or r m ) to be "remembered” for some distance downstream; the 
distance is dependent on the particular flow situation. Experience has shown that two 
stream flows with u a /uj < 0.25 can be started with a "fully developed" shear stress. 
This fully developed shear stress is obtained from the corresponding fully developed 
two-dimensional shear layer. Even though the present method is generally limited to 
flows with thin initial boundary layers, two stream jet flows with u a /uj < 0.25 and very 
thick initial boundary layers (jet nozzle boundary -layer thickness up to 0.4r o ) have been 
successfully computed with the following procedure: (1) the initial boundary layer is 
assumed to be negligible, and (2) the inner shear layer radius q is adjusted to match 
the experimental value of the excess momentum. 

Even for thin boundary layers, the influence of the initial conditions persists 
throughout the flow field when u a /Uj exceeds about 0.3. Therefore, the concept of a 
negligible initial boundary layer and a fully developed initial shear stress is not usable 
for such flows. In order to treat these flows, a control volume analysis of the initial 
region has been developed. 

Control Volume Analysis of Initial Region 

A sketch of the initial region just downstream of the nozzle lip is shown in figure 8. 
The initial boundary layers are characterized by power law velocity profiles: 

a 


and 

u_ ( I o - r 

U 3 V 6 i 

The nozzle wall is assumed to be adiabatic; therefore, the stagnation temperature is 
constant in each boundary layer. Specification of the wall skin-friction coefficients, 
c fa and c^, and a^ = 0.3 completely defines the mean flow and the turbulence quanti- 
ties at the initial station b£. The wall skin-friction coefficients are determined with the 
method of Spalding and Chi (ref. 21). 
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At some downstream station, si, the flow is assumed to have evolved to the fully 
developed shear layer profile shapes for velocity and near -field TKE. The following 
assumptions are made about the process between bi and si: 

(1) There is no net entrainment into the shear layer. 

(2) The process occurs at constant pressure. 

(3) The excess momentum is conserved. 

(4) The length scale of the process is sufficiently small to insure that the volume 
integral of turbulent dissipation is negligible. 

The following equations are written between stations W and si: 

Momentum: 


f* r o + ^a 

p r i+b 

\ pu(u - u a )r dr 

J 0 

= J 0 PH U ~ u a) r dr 


Continuity: 


(15) 


f r o+6 a 

pri+b 

\ pur dr 

= \ pur dr 

- 5* 

0 °J 

U r i 


Turbulent kinetic energy: 


(16) 



With 6 a , 6j, n a , nj, Cj a , and Cjj specified, equations (15) to (17) are solved for rj, 
b, and k m . These values are then used as initial conditions for the integral TKE solu- 
tion of the remainder of the flow field. 


This control volume analysis is obviously not applicable as u a /Uj approaches 
unity, since in such flows a very large distance is required to approach a fully developed 
profile shape. In addition, the initial region analysis as formulated is not applicable to 
mixing of streams with greatly different densities. The computations are made with 
K = 1 at station si, and the solution does not, in general, conserve species or energy. 
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RESULTS AND DISCUSSION 

Initial conditions for the computed experimental flows are presented in table I. 

The initial shear stress levels are characterized by the turbulent Reynolds number R-p. 
In addition to being used to specify the initial shear stress, the axial distribution of R-p 
throughout the flow field is very informative, and a number of such distributions is pre- 
sented. For the cosine profile, the midpoint shear stress in the mixing layer is given by 



p m (Au) 2 2R T 


where Au is the local velocity difference across the layer. If self preservation is 
approached, then R T must become constant in the flow field. 

Most of the flows considered in this conference have a constant -pressure boundary 
condition, whereas the analysis was developed for a ducted system with a prescribed duct 
wall shape. The constant -pressure axisymmetric flows were computed in a very large 
cylindrical duct (r w = 1000r o ) so that negligible axial pressure gradients were predicted. 
Integrated momentum in the duct is conserved to a high degree of accuracy in the com- 
putations, typically to within one part in 10 5 . The degree to which the jet excess momen- 
tum is conserved in constant -pressure flows is illustrated by the fully developed axisym- 
metric jet (test case 18). For similar velocity profiles, the product of mixing zone width 
and center-line velocity should remain constant throughout the second regime. This 
product changed 1.5 percent from the value at the end of the core at x/D = 50, and 
2.8 percent at x/D = 100. Even though the excess jet momentum is not exactly conserved 
in the calculations because of the ducted boundary condition, the results are considered to 
conserve excess momentum adequately when the inprecision of most experiments is taken 
into account. 

Effect of Velocity Ratio on Growth of Fully Developed Two-Dimensional 

Shear Layer - Test Case 1 

As presently formulated, the integral analysis cannot be used for zero secondary 
velocity. Therefore, all computations of flows with nominal zero secondary velocity 
were made with a secondary velocity 0.01 times the maximum velocity in the flow field 
(Uj or u Q ). The fully developed two-dimensional shear layers were computed as thin 
shear layers at the boundary of a large axisymmetric jet; the duct radius r w was set 
equal to 100r o for these cases. In no case did the predicted shear-layer thickness 
exceed 2 percent of the jet radius. 

Computations for test case 1 were made for u 2 /Uj = 0,01, 0.2, 0.4, 0.6, and 0.8. 

The results, shown in figure 9 (test case 1), fall on the classic relationship 
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°o _ U 1 u 2 
a U j + u 2 


where a Q = 12.9 for u 2 /uj = 0. The standard relation for a used in this computa- 
tion becomes 


a 


3.12 

db/dx 


for the cosine profile. It should be noted that a = 12.9 by this definition corresponds 
to o=12 when the cosine profile midpoint slope is matched to the widely used error 
function profile. The fully developed R T varied negligibly from 143 over the entire 
range of u 2/ /uj. 


Effect of Mach Number on Growth of Fully Developed Two-Dimensional 

Shear Layer — Test Case 2 

Results for test case 2 are shown in figure 10(a). The computations were made with 
u 2/ u l = therefore, all o values are about 2 percent too large. The ratio <r 0 /cr, 
where a 0 is the value at Mj = 0, is shown in figure 10(b), along with the value of fully 
developed R T . Based on available experimental information on a, the predicted a Q /a 
probably decreases too abruptly in the Mi range of 0.5 to 1.5. The predicted a values 
in this Mi range can be altered by slight refinements of the a 2 function in the appro- 
priate Rt range. The predicted a values are considered to be good at Mi = 2 and 
Mi = 3. 


Effect of Density Ratio on Growth of Fully Developed Two-Dimensional 

Shear Layer - Test Case 3 

The computations for test case 3 were made for low-speed flow and u 2 /Uj = 0.2. 
Results for cr and R-p are shown in figure 11 as a function of P\/P2‘ Evaluation of 
the cr predictions at high Pj/P 2 is nearly impossible because of the lack of experi- 
ments in this range. The cr results for this case are not influenced by the factor which 
causes the density ratio, that is, temperature difference or molecular weight difference. 

In general, very large axial distances were required to approach the fully devel- 
oped condition in all these shear -layer computations. All the flows were computed for 
an axial distance of several hundred initial shear-layer thicknesses; such distances were 
required to approach closely the fully developed condition unless the initial shear stress 
was luckily chosen to be very close to the fully developed value. These results clarify 
the extreme difficulty in accomplishing a shear-layer experiment in which the flow truly 
approaches a fully developed condition. 
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Maestrello and McDaid Axisymmetric Jet - Test Case 6 

The computations for test case 6 were made in two ways. In the first (curve a, 
fig. 12), the experimental profile at x/r Q = 2 was fitted with a cosine profile, and the 
computations were started with a fully developed shear stress. The predicted rate of 
mixing is too large; therefore, the shear stress does not yet approach the fully developed 
value at x/r Q = 2 but is somewhat lower. The computations were also started at x = 0, 
with negligible initial boundary layer, fully developed R T and r i corrected to achieve 
the excess momentum shown at x = 2r 0 „ This second computation (curve b, fig. 12) 
yields a first regime which is somewhat too long, but the results are better than those of 
the first computation. 

The abrupt change in center-line velocity at the end of the first regime is charac- 
teristic of the integral method but is of little concern unless the main interest is in the 
transition region at the end of the potential core. 

Eggers Supersonic Jet Into Still Air - Test Case 7 

The prediction of this well-defined experimental flow (test case 7) is very satis- 
factory. (See fig. 13(a).) Computations were started by assuming a negligible initial 
boundary layer and fully developed R T (from fig. 10(b)). The predicted potential core 
is slightly longer than that shown by the experiment, but the far -field agreement is excel- 
lent. The predicted velocity profiles at x/r Q = 8, 27, and 99 are also satisfactory. 

(See fig. 13(b).) These profiles illustrate that the cosine profile approaches zero at the 
outer edge of the layer more rapidly than does the experiment; the cosine profile is gen- 
erally better for two stream flows. 

The predicted axial variation of Rt (fig. 13(c)) shows a very large change in R T 
for this flow (from 283 to 82). Thus, most of the a 2 - Rt function (fig. 7) was used in 
this prediction. 


G.E. Heated Subsonic Jet - Test Case 8 

The experimental velocity profile at x/D = 2.79 was fitted with a cosine profile; 
the computations were started with this velocity profile and fully developed Rt- Pre- 
dictions of both center -line velocity and center -line static temperature are very satis- 
factory for test case 8. (See fig. 14.) 

For stall Jet in Moving Stream - Test Case 9 

This flow (test case 9) was computed by assuming a thin initial boundary layer and 
fully developed Rt; r i was corrected to yield the experimental excess momentum and 
the actual duct radius r w /r 0 = 16 was used for the prediction. The predicted center- 
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line velocity agrees very well with the experiment at all axial locations. (See fig. 15(a).) 
The detailed far -field behavior of the excess center-line velocity is illustrated in the log 
plot (fig. 15(b)). Because of the assumption of unity Prandtl and Lewis numbers, the pre- 
dicted a is identical to the excess velocity (u c - u e y(u 0 - u e ). Predicted half- 
velocity width (fig. 15(c)) agrees well with the experiment for x/D up to 25, but falls 
about 10 percent under the experiment at x/D = 80. 

Chriss Hydrogen Jet in Moving Air Stream - Test Case 10 

This flow (test case 10) was computed in two ways. In the first, a cosine velocity 
profile was fitted to the data at x/D = 2.97. This profile, along with the experimental 
shear stress shown in reference 22 (R x = 98.3), was used to start the computations. 
Predicted results for the center-line velocity are fairly good (fig. 16(a)), but the pre- 
dicted rate of decay of center-line concentration is too low. This type of concentration 
prediction is typical of the integral analysis, since it is limited to unity turbulent Prandtl 
and Lewis numbers. 

The second computation was started at x = 0 with small initial shear layer thick- 
ness and fully developed shear stress (R t = 120). Predicted center -line velocities 
(fig. 16(b)) are somewhat better than those for the first computation, but, of course, the 
center -line concentration decay rate is again underpredicted. 

Eggers and Torrence Axisymmetric Jet in 
Moving Air Stream - Test Case 1 1 

This flow (test case 11) has wakelike behavior and it is unlikely that it can be prop- 
erly predicted by a constant -pressure mixing analysis. In spite of the thick initial bound- 
ary layers and the velocity ratio (u e yiij = 1-36), computations were started at x = 0 with 
thin initial shear layer and fully developed shear stress. The predicted center-line 
velocity distribution is surprisingly close to the experiment. (See fig. 17.) Wakelike jet 
flows, with jet momentum flux less than external stream momentum flux, tend to have 
relatively high relative shear levels (low R T ), and it is possible in such flows that the 
influence of the initial conditions does not persist very far downstream in the flow field. 
This aspect should certainly be further explored. 

Eggers Hydrogen Jet in Moving Air Stream - Test Case 12 

As with the preceding case, this flow (test case 12) has jet momentum flux less than 
external stream momentum flux. In this case, however, the jet momentum deficit is 
caused by density rather than by velocity ( u e/ u j = 0.37j. Again, the computations were 
started with a thin initial shear layer and fully developed R T . Although the length of the 
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first regime is somewhat overpredicted (fig. 18) , the downstream prediction for both 
center -line velocity and center-line concentration is very good. In view of the unity 
Schmidt number assumption, the good prediction of both velocity and concentration was 
unexpected; the results are undoubtedly caused in part by the neglect of the initial bound- 
ary layers. 


Chevray Axisymmetric Wake — Test Case 15 

Computation of this flow (test case 15) was started at x = 0 with an equivalent 
second regime cosine velocity profile; the cosine profile was selected to match the experi- 
mental momentum defect and mass flow in the shear layer. An initial value for R T was 
established by equating the experimental TKE flux (ref. 23) with the TKE flux in the equiv- 
alent fully developed profile. The theory does not predict the initial rapid acceleration of 
the center-line velocity (fig. 19(a)) for x/D < 2. This rapid acceleration is caused in 
part by a favorable pressure gradient just downstream of the body; the axial pressure 
gradient was neglected in the calculations. For x/D > 4, the center-line velocity defect 
W is underpredicted; however, the log plot of W against x (fig. 19(b)) shows that the 
predicted decay rate at x/D = 18 is somewhat less than the experimental rate. The 
predicted axial distribution of W follows an x _1 decay from x/D = 5 to the maximum 
axial distance computed (x/D = 200). Better prediction of the decay rate at x/D = 18 
would require less dissipation than that used in the calculations. On the other hand, 
achievement of the x -3 / 3 decay rate for W, as predicted for self-preservation, would 
require much higher dissipation. One can only conclude that (1) the far -field dissipative 
mechanism is much different from the near -field mechanism, or (2) self -preserving axi- 
symmetric wakes are never attained. The lack of far -field experiments on the axisym- 
metric wake makes it difficult to decide which of these conclusions is correct. 

The predicted axial variation of R-j. (fig. 19(c)) shows that Rx never approaches 
a constant value but continuously decreases in the axial direction. 

Demetriades Supersonic Axisymmetric Wake — Test Case 17 

This flow (test case 17) was computed by fitting a cosine profile to the experimental 
velocity data at x/D = 17, and by using the Rq> value quoted by Demetriades (ref. 24) 
for x/D * 17. The plot of — against x/D (fig. 20(a)) shows that the experimental 

W 3 /2 

center -line velocity defect is well predicted. A log plot (fig. 20(b)) shows that the pre- 
dicted W decays even more rapidly than an x"* decay. Again, this axisymmetric 
wake prediction is very different from the x -3 / 3 decay of W which is required for 
self-preservation. The predicted axial variation of Rx (fig. 20(c)) is considerably dif- 
ferent from that shown by Demetriades (ref. 24); he showed Rx to be nearly constant 
at 32 for x/D greater than about 30. 
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Fully Developed Axisymmetric Jet - Test Case 18 

This important and fundamental flow (test case 18) was computed from x = 0 by 
assuming a thin initial shear layer and a fully developed near -field shear stress 
(Rt = 143). The predicted center-line velocity (fig. 21(a)) follows an x -1 decay for 
x/D greater than 15. Although the predicted center-line velocity agrees well with the 
experiment of Albertson et al. (ref. 25) at all axial stations, the agreement with the 
Wygnanski and Fiedler experiment is fairly satisfactory only for x/D > 40. The 
Wygnanski and Fiedler flow apparently does not approach self-preservation as quickly 
as do other reported jet experiments; the terms "fully developed" or "self -preserving" 
are questionable when applied to the Wygnanski and Fiedler experiment. 

The predicted velocity profile (fig. 21(b)) is based on the local mixing-zone growth 
rate at x/D = 60. Compared with the experimental profile of Wygnanski and Fiedler, 
the half -velocity radius is well predicted, but the predicted profile is fuller near the 
center line; the predicted profile is closer to the empirical profile of Albertson et al. 

(ref. 25) in the high-speed half of the shear layer. The theoretical profile approaches 
zero velocity in the outer part of the shear layer more rapidly than does the experi- 
mental profile. The predicted TKE level (fig. 21(c)) agrees well with the Wygnanski and 
Fiedler experiment in the outer half of the profile, but the predicted center -line TKE 
level is about 13 percent below the experiment. It should be noted that the predicted TKE 
level is satisfactory in the region of peak shear stress. Predicted axial variations of 
Rt and the TKE profile parameter w are shown in figure 21(d). Even though the 
center-line velocity closely follows the x~l decay (required for self-preservation) 
downstream of x/D = 15, R'j and w do not become constant until considerably far- 
ther downstream. 


G.E. Heated Supersonic Jet - Test Case 19 

This flow (test case 19) was computed by fitting a cosine profile to the velocity data 
at x/D = 2.79, and by using a fully developed initial shear stress. Agreement between 
predicted and experimental center-line velocity (fig. 22) is fairly satisfactory. 

Paulk Jet in Moving Stream - Test Case 20 

Even though the initial boundary layers are fairly thin in this flow (test case 20), the 
velocity ratio (ue/u 0 = 0.48) exceeds that where a fully developed initial shear stress can 
be used. A potential core length of approximately 12D is predicted with the assumption of 
a fully developed initial Rt- The control volume initial region analysis was applied with 
6 a = 0.082r o and 6j = 0.063r o . The experimental boundary -layer shapes are not well 
defined in reference 22; therefore, the initial region computations were carried out for 
two values of the velocity profile exponent (n a = nj = 1/4 and n a = nj = 1/7). The result- 
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ing values of b, r^, and shear stress were then used to start (at x - 0) the computa- 
tions for the entire flow field. Both computations are shown in figure 23; the predicted 
flow field is obviously not very sensitive to the initial boundary -layer shape. Both pre- 
dictions agree fairly well with the experimental center -line velocity but, of course, the 
experimental center -line concentration decays more rapidly than predicted. 

It should be noted that the experimental values of (u c - u e )/(u 0 - u e ) were com- 
puted with u 0 = 390 ft/sec (119 m/sec) and u e = 187 ft/sec (57 m/sec), the observed 
values of u 0 and u e downstream in the flow field. 

Chriss Hydrogen Jet in Moving Air Stream - Test Case 21 

As with test case 10, this flow (test case 21) was computed (1) with experimental 
initial conditions , and (2) with fully developed R-j- at x = 0. The initial shear stress 
for the first computation was taken from reference 22 at the most upstream axial station 
for which a shear stress was measured in this flow. The predicted center -line velocity 
agrees well with experiment (fig. 24(a)) but, as usual, the center-line composition decay 
is underpredicted. The computations which were started at x = 0 are shown in fig- 
ure 24(b); the agreement with experiment is poor. Therefore, it can be concluded that the 
velocity ratio of this flow (u e /u 0 = 0.31) exceeds that for successful use of a fully devel- 
oped initial shear stress. 

Eggers Hydrogen Jet in Moving Air Stream - Test Case 22 

As with test cases 11 and 12, this flow (test case 22) is wakelike in that the jet 
momentum flux is less than the outer stream momentum flux. The flow also has very 
thick initial boundary layers. In spite of the thick boundary layers and the velocity ratio 
( u e/ u j = 0.54), computations were started at x = 0 with a thin initial shear layer and 
with fully developed R'p . The predictions (fig. 25) are not as bad as could be expected, 
but the analysis does not predict qualitatively the behavior of the center -line velocity. 

It is unlikely that this flow will ever be properly predicted with an analysis which ignores 
the significant pressure gradients in the flow field. a 

CONCLUDING REMARKS 

When the rather small investment of empirical information on turbulence structure 
is considered, the quality of the mean flow predictions is surprisingly good. Experience 
in computing a variety of flows has not indicated that the assumption of universal TKE 
profile shapes is grossly incorrect. The utility of the admittedly oversimplified model 
for turbulent dissipation has been significantly increased by relating the dissipation coef- 
ficient a 2 to the local turbulent Reynolds number. 9 * 
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A persistent point of criticism about TKE methods for free shear flows has been 
that the required detailed information on initial conditions is not generally available for 
engineering flows. Experience has shown that a considerable number of practical flows 
can be computed without detailed information on the initial conditions. For u e /u Q less 
than 0.25, specifying a fully developed initial shear stress works very well. Wakelike 
flows, such as test case 12, seem to be fairly insensitive to initial conditions. For jet 
flows with velocity ratio u e /u 0 greater than 0.25, the control volume initial region 
analysis shows promise. The very limited experience with this initial region analysis 
for constant density flows has indicated that the results are not overly sensitive to the 
initial boundary -layer characteristics; good guesses about the initial boundary layers 
seem to be sufficient. 

Turbulent kinetic energy methods, of which this method is a simplified example, 
represent a fundamental improvement over eddy -viscosity models in that more of the 
physics of the turbulent motion is taken into account, albeit crudely. It appears that 
these methods have been developed to the point where they are routinely applicable to 
engineering calculations on a broad class of free turbulent flows. 
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TABLE I.- INITIAL CONDITIONS FOR THE EXPERIMENTAL FLOWS 


Test 

case 

Flow 

X 1 

b l 

r u 

r T,1 

Notes 

6 

Maestrello and McDaid axisymmetric jet, 

2r 0 

0.87r o 

0.52r o 

190 

Experimental profile, fully developed Rx 


Mq = 0.64 

. 0 

•01r o 

,86r 0 

190 

r Q corrected for momentum 
Fully developed R T 

7 

Eggers supersonic jet, M 0 = 2.22 

0 

•01 r„ 

r o 

283 

Fully developed R T 

8 

G.E. heated jet, M 0 = 0.7 

5,58r 0 

1.70r o 

.18r 0 

212 

Experimental profile, fully developed R T 

9 

For stall jet in moving stream, 
u e /u 0 = 0.25 

0 

•01 r o 

o 

u 

CO 

CD 

143 

r 0 corrected for momentum 
Fully developed Rx 

10 

Chriss H 2 jet in air stream , 

5.94r 0 

1.132r 0 

M2t q 

98.3 

Experimental profile and Rx 


u e/ u o = °- 16 

0 

•01r o 

T o 

120 

Fully developed R T 

11 

Eggers and Torrence jet in moving stream 

0 

.0lr o 

r o 

130 

Fully developed R T 

12 

Eggers H 2 jet in air stream, M e = 1.33 

0 

•01r o 

r o 

125 

Fully developed Rx 

15 

Chevray axisymmetric wake 

0 

.0257D 

0 

1930 

Equivalent cosine profile, u c = 0.083u e 

17 

Demetriades axisymmetric wake, 
M e = 3 

17D 

1.00D 

0 

48 

Experimental profile. Rx quoted 
by Demetriades 

18 

Fully developed axisymmetric jet 

0 

.01r o 

r o 

143 

Fully developed R T 

19 

G.E. heated jet, Mq = 1.36 

5.58r 0 

1.30r o 

•46r 0 

238 

Experimental profile, fully developed Rx 

20 

Paulk jet in moving stream, 

0 

.113r 0 

.916r 0 

36.8 

Initial condition analysis, n a = n^ = 1/4 


u e /u 0 = 0.48 

0 

.122r 0 

,923r 0 

42.9 

Initial condition analysis , n a = nj = 1/7 

21 

Chriss H 2 jet in air stream, 

5.15r 0 

1.17r 0 

.464r 0 

48.2 

Experimental profile and Rx 


u e /u 0 = 0.31 

0 

•01r o 

r o 

130 

Fully developed Rx 

22 

Eggers H 2 jet in air stream, M e = 2.5 

0 

.01r o 

r o 

110 

Fully developed Rj 
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Figure 5.- Family of turbulent kinetic energy profiles. 
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Figure 6.- Center-line turbulent kinetic energy for axisymmetric jet into still air. 
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Figure 7.- Dissipation coefficient as a function of turbulent Reynolds number. 
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Figure 8.- Nomenclature for initial region analysis, si and hi denote 

shear layer and boundary layer. 
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Figure 9.- Test case 1 - Effect of velocity ratio on growth of constant-density 
two-dimensional fully developed shear layer. 
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Figure 10.- Test case 2 — Effect of Mach number on fully developed spreading parameter 
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Figure 12.- Test case 6 — Maestrello and McDaid axisymmetric jet into still air 
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Figure 22.- Test case 19 - Heated axisymmetric jet into still 


air. M 0 = 1.36. 
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Figure 23.- Test case 20 - Paulk axisymmetric jet in moving stream. u e /u 0 = 0.48. 
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Figure 24.- Test case 21 - Cfyriss H 2 jet in moving stream. u e /u 0 = 0.31. 
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DISCUSSION 


S. W. Zelazny: You assume that the Schmidt number was unity in your analysis. Wouldn’t 
it be a simple exercise to remove that restriction by introducing an additional equation? 

C.E. Peters: Yes sir, we plan to do that. It’s a Uttle frustrating because obviously I 
have access to Tom Harsha, who sits 15 feet away, and for those flows where I’m inter- 
ested in a nonfully developed case or nonunity Schmidt number I can always use the finite- 
difference results. But there are certain flows that I can handle with the integral method 
we can’t handle yet with the finite -difference method. So I will probably go ahead and try 
to do something to approximately handle nonunity Schmidt number. 

S. C. Lee: You mentioned this method you have has some advantage over finite-difference 
schemes. 

C.E. Peters: I didn't refer to it as an advantage, I called it a difference. It's a different 
set of approximations. Instead of assuming an aj function we just assume that the tur- 
bulent kinetic energy profile shapes are invariant as a different approach. 

S. C. Lee: Perhaps you could summarize for us what would be the difference between 
this method and the finite-difference method as Harsha presented it. 

C.E. Peters: The point is that my kinetic-energy profiles are universal in the given 
position in terms x/x core . Because of the variable -density and variable -velocity ratio, 
the shear profile varies from flow to flow; that is, the normalized shape of the shear pro- 
file. Therefore, a^ varies in shape from flow to flow, a condition implied through the 

whole analysis. I backed it out for some flows, and it does vary from flow to flow. So 
it’s an alternate approximation to saying something about aj. That's the only fundamen- 
tal difference. Tom and I are using essentially the same dissipation, and, of course, in 
return for my specifying the profile shape of kinetic energy, I don't have to say anything 
about diffusion. 

M. V. Morkovin: I think this is another case where we can benefit from the comparisons 
of different methods. Could you tell us something about your general experience, when 
you prefer to use Harsha' s method to yours and how much more complex it is. I gather 
you do have some problems in sensitivity to initial conditions. If you could give us a 
little bit of briefing on that, I think it would be highly valuable. 

C.E. Peters: I'll give you my experience first and Tom can comment on it if he likes. 

The implication of my being able to use a fully developed initial shear stress up to a 
velocity ratio of 0.25 merely says that these flows are strong shear flows, and they wipe 
out any small inaccuracies in the initial condition fairly quickly. So that seems to be a 
viable procedure. When the flows become wake-iike in nature, of course, the initial con- 
ditions persist indefinitely; in that case my analysis is not applicable at all. For example, 
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if the hole never fills up, it's not useful. So I limit it to a velocity ratio of maybe 0.7 or 
0.8 as the upper limit. I haven’t tried to define carefully an upper limit here. Certainly, 
the finite-difference method is much more powerful in a sense of handling nonsimilar ini- 
tial conditions and so forth. It's interesting to note the sort of generalized information 

one can back out with an alternate set of assumptions. For example, about a x profiles. 
So in addition to the practical requirement that we have a better transport model for our 
"nuts and bolts” engineering calculation schemes, we are also interested in some of the 
implications in general. 

Anon. ; Which costs more? 

C^E. Peters: This is embarrassing for me because Harsha's program runs faster. But 
one of the reasons is that this program is a mess; it does lots of different things. This 
regime is far beyond what we were talking about in this particular calculation. We even 
do circulating flows with them by using patched up techniques. So the program is not 
optimized for this kind of flow. If we build a constant- pressure version, I am sure we 
can get our run time down perhaps an order of magnitude. 

— — Alber : Your calculation of the spreading parameter variation with Mach number for 
the two-dimensional shear layer shows quite an increase in the variation of sigma with 
M. That is your figure 10. Now in this same calculation by Harsha, which I believe is 
a very similar model to yours, he shows hardly any variation at all. 

C^E. Peters: We can’t explain that - it’s different. Tom, would you care to comment? 

— T ~ H Hghg: : Well, the only comment I could make really is that with shear layers I was 
forced to use an almost patently ridiculous &1 model which Peters does not have to use 
because his a^ comes out of shear stress profile and kinetic energy profiles. I would 

suspect that the problem is that the model that I used is simply inappropriate for the 
shear layer. 

J. M. Eggers: I was wondering if you or anyone else here could comment on what we 
could do to improve the modeling in the transition region to get rid of this somewhat 
atrocious inflection point at the end of the core region. 

C. E. Peters : The sharp corner on the center-line velocity distribution is characteristic 
of integral methods, at least as we have put them together. It doesn’t bother me very 
much because most of the required information is not centered in that particular region, 
but the finite difference smooths it out. That’s the idea of patching together two regimes 
dis c ontinuous ly . It’s bound to give a sharp corner with an integral method without some 
relaxation of profile shapes locally or something. 

— B - Spalding: I have a question and a comment about R T . The question is a simple 
one - in the paper R T seems to be defined twice, once it has the difference between 
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u c and u a on the top and in the second it has the square of those. Which of those was 
actually used ? 

C. E. Peters : I think the difference comes because when you plug in the midpoint velocity 
profile slopes, it is in terms of delta u over width and I think that's where it comes from. 

D, B. Spalding : The other remark I want to make is that this quantity appears to be 

related to the ratio of the production term to the dissipation term which Rodi found also 
to be significant, and it will be interesting when we have time just to compare whether 
there was a quantitative connection between your function of and our function of 

production over dissipation. 

C. E. Peters: Yes, I certainly agree. We haven't had time to give much thought as to 
what this R T function means in terms of structure, but it would be nice if we could at 
least rationalize, in terms of energy spectrum or whatever, why the length scale changes 
or the R t changes the effective dissipation. 

D. B. Spalding : At least you can see that with your definition, when u c equals u a , 
there will be no shear stress, so production is zero. And so R-j- is zero; it’s the same 
as Rpy/e is equal to zero. Even closer quantitative connections are being worked out. 

C. E. Peters: It was also interesting for me to read in Joe Schetz's written version that 
the constant in his eddy-viscosity model perhaps should be related to a similar parameter 
to this - the ratio of U prime over delta U. So, this parameter is obviously developed 
by mean flow correlations during the near field of the wire and jet series, and I just 
stumbled across it. I think it is better than a mean flow parameter, such as Mach num- 
ber or some density parameter, and I think it is related to the turbulent structure. 

S. Corrsin : I was interested in R T because the numbers look so much larger than the 
ones I computed 15 years ago.* Did you compute it for a wake also? 

C. E. Peters: Yes, it was shown in the paper I think on the Chevray case. Remember my 
length scale is the full width of the shear layer; that is, from the centerline out to the 
outer edge. 

S. Corrsin : Well, I was just looking at my paper from 1957 in which is used the half- 
width and the full- velocity difference. I got a value of 12 for Townsend's wake. 

C. E. Peters : I didn't do the two-dimensional wake. 

S. Corrsin : And for the round jet, I defined it in terms of the momentum diameter and 
it came out to be 15. 


* Corrsin, S.: Some Current Problems in Turbulent Shear Flows. Symposium on Naval 
Hydrodynamics (Washington, D.C.), Sept. 24-28, 1956, pp. 373-400. 
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C.E. Peters : My number of 75 for the asymptotic value for the fully developed jet is con- 
sistent, considering the difference in length scale definition, with the number tabulated 
experimentally by B. G. Newman. 

S. Corrsin : By whom ? 

C. E. Peters ; Newman, in his survey paper a few years ago. So I think it is consistent 
with other quoted information. 

S. Corrsin: I can make a physical comment about this concept. I first heard it suggested 
back in the early 1940's by Hans Liepmann, who proposed the idea that perhaps turbulent 
shear flows tend to keep themselves in a state of lower critical Reynolds number based 
on the turbulent viscosity, and it's a sort of self-destroying system that always disturbs 
itself violently. That was sort of the reason that I computed these. For bounded flows, 
they seem to vary more, but the general idea was that for shear flows without boundaries 
there are probably universal constants for each geometry. 
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